Stationary and transient leakage current in the Pauli spin blockade 
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We study the effects of cotunneling and a non-uniform Zeeman splitting on the stationary and 
transient leakage current through a double quantum dot in the Pauli spin blockade regime. We find 
that the stationary current due to cotunneling vanishes at low temperature and large applied mag- 
netic field, allowing for the dynamical (rapid) preparation of a pure spin ground state, even at large 
voltage bias. Additionally, we analyze current that flows between blocking events, characterized, in 
general, by a fractional effective charge e*. This charge can be used as a sensitive probe of spin 
relaxation mechanisms and can be used to determine the visibility of Rabi oscillations. 
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Initialization and readout of well-defined quantum 
states are necessary for spin coherence measurements [1], 
single-spin resonance [2, 3, 4], and quantum information 
processing. Single electron spins in quantum dots show 
promise for quantum information tasks [5] due to their 
long coherence times, but their quantum states can be 
difficult to initialize (relying on slow spin relaxation pro- 
cesses) and read out. The Pauli spin blockade (PSB) [6] 
partially solves these problems, where current through a 
double quantum dot (DQD) is blocked conditional on the 
microscopic spin state of electrons. 

The PSB is, however, imperfect; hypcrfinc interaction 
between electron and nuclear spins in III-V semiconduc- 
tors lifts spin selection rules and can lead to a finite 
leakage current [7, 8]. Very recently, PSB has been ob- 
served in DQDs made from silicon [9, 10] and carbon 
nanotubcs [11, 12], in which the majority isotope has no 
nuclear spin. Even in these systems, the PSB can be 
lifted through spin exchange with the leads due to, e.g., 
cotunneling processes [13, 14, 15]. Significantly, cotun- 
neling events have been shown to be essential even in 
nuclcar-spin-carrying quantum dots to describe nuclear- 
spin polarization in the PSB regime [16, 17], and there- 
fore should be taken into account. 

Single-spin resonance measurements often rely on the 
transient current that flows before current is blocked as a 
probe of the electron spin state [2, 3, 4]. An anomalously 
large transient current has recently been found [2] , char- 
acterized by an effective charge e* that passes through 
a DQD between blocking events. Without a complete 
understanding of this additional leakage, it may not be 
possible to determine the visibility of Rabi oscillations in 
these systems. 

Here, we evaluate both the stationary and transient 
leakage current through a DQD, giving simple analytical 
expressions for the stationary current and the transient 
effective charge e*. We find that e* reaches universal 
fractional values and that a measurement of e* in general 
can be used to extract valuable information related to 
slow spin relaxation processes. 
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FIG. 1: (Color online) A DQD coupled to leads (a). The 
charge stability diagram is shown in (b). At low tempera- 
ture, inelastic cotunneling (c) induces transitions to lower- 
energy dot levels. Two-electron states at e = are shown 
in (d) with allowed T = Q cotunneling transitions indicated 
with straight dashed arrows and curved arrows indicating 
sequential-tunneling processes. 



We consider a series-coupled DQD in a magnetic field 
gradient (Fig. 1(a)). A field gradient is important in 
spin resonance experiments for local addressing, and may 
arise from the stray field of a nanomagnet [4] or the Over- 
hauser field due to non-uniformly polarized nuclear spins 
[2]. We work in a regime where only the (0, 1), (1, 1), and 
(0, 2) charge states are relevant (the triangular region in 
Fig. 1(b)). Here, {riL.nB) indicates m electrons in dot 
orbitals / = The Hamiltonian in this projected 

subspace is H ^ Hi^ + V, where 

= Y. 1^') I + E ^lkcl.^lka , (1) 
j Ika 

V = 5b\S){To\ + tiAllcika\j){3'\+\^.c. (2) 

klajj' 

The first term in iJo describes the eigenstates |j) = 
{[(t), |a)} of the unperturbed DQD, with single- electron 
states labeled by spin a = {Tj i} and two-electron states 
labeled with a, shown in Fig. 1(d) (three spin triplets. 



(a = r„ 



0,±) and two spin singlets, (a = S±), 



giving hybridized (1,1) and (0,2) charge states due to 
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FIG. 2: (Color online) Leakage current in the PSB regime 
due to inelastic cotunneling processes. We have taken e = 0, 
Fi = Tfl = lO/iel/, t = 100 ^tel/, A = Imel/, and g = 2.0. 
We show the evolution of /(-B) as &B is varied from mT 
(dashed line), to 20 mT (dotted line), 200 mT (solid line), 
and for Wtq^o- S> Wa^f3, T = (dash-dotted line, from Eq. 
(7)). Evolution of the curve for SB = 200 mT as temperature 
is raised is shown in the inset. 



an interdot tunnel coupling t [18]). The second term in 
Hq gives the energy of Fermi-liquid leads. The Zeeman 
gradient 6b = g^B{BL — Br)/2 couples the (1, l)-singlct 
15*) and ttt-s — triplet \Tq) (here, g is the g-factor and Bi 
is a local magnetic field in dot /). The second term in V 
describes hopping processes from dot I to lead I with cou- 
pling ti and matrix elements = {i\d\^\j'). Here, 
creates an electron in dot I with spin a and c\j^^ creates 
an electron in lead / and orbital fc, with spin a. 

Working from a standard Hubbard model for the DQD, 
we find the energies -E^kx) = — (e A)/2 -f- (— )&/2, where 
e is the detuning (energy difference) between the (1,1) 
and (0, 2) charge states and A controls the depth of the 
two-electron levels (see Fig. 1(c)). For the two-electron 
states, we have Et^ = —A, £'t± = —A ± fe, and E'g^ = 
-A - e/2 ± Ve^ -1-8*2/2, with h = gusiBL + Bb)/2. 

We solve for the diagonal elements of the reduced 
(DQD) density matrix pi with the Pauli master equa- 
tion: 

P^ = Y.^P,W,^^-P^W,^J)■ (3) 
3 

We calculate the transition rates Wi-,j = Wf\,j + Wf^^ 
directly from Fermi's golden rule. Here, the sequential- 
tunneling rates {W^^^ = Y^iW^^^ oc |*;p) describe di- 
rect hopping at leading order in the dot-lead coupling and 
the cotunneling rates (W™*^- oc \ti\^ ^ rj > 2) are higher- 
order in ti. 

We consider standard initial conditions, with lead / 
held in thermal equilibrium with Hamiltonian i?o at 
chemical potential At first order in y, only the sec- 
ond term in Eq. (2) contributes to transport, giving the 



usual sequential-tunneling rates (fi = 1) [19]: 

a' 

= Y.^i\AZ\^\i-hM\. (4) 

a' 

Here, F; = 2'Kv\ti\^ ^ where v — 'Ylik^^^P ^ ^'fc) 
density of states per spin at the Fermi level, /; {E) is a 
Fermi function at temperature T and chemical potential 
[ii, and = - Ej. 

For large bias, /i^ — /i^ = A/.J > |A| > 2kBT, the 
stationary current is given by / = e J^a ^a-^aPa, where 
Pi is a solution to Eq. (3) with pi = 0. At leading order 
in V, current will be blocked if one of the triplet states is 
populated, since = 0. This is the PSB effect. In 

the absence of other spin-relaxation mechanisms, higher- 
order contributions in V must be considered to explain a 
finite leakage current. At second order in V wc find: 

W^To^. = [1 - fBi^To.)] , (5) 

9F2 // 

^ ' aa'a" 

where F(a;,r) = w/(l-e-"/'=«'^). Eq. (5) gives the rate 
for field-assisted sequential-tunneling processes, where Sb 
converts To to a singlet, which can then escape from the 
DQD via first-order (sequential) tunneling to the right 
lead (see the curved dashed arrow in Fig. 1(d)). Eq. 
(6) gives the rate for an inelastic cotunneling process 
(Fig. 1(c)), allowing for conversion of triplets to singlets 
(dashed straight arrows in Fig. 1(d)). A competition be- 
tween the two rates in Eqs. (5) and (6) will determine 
the leakage current in the PSB regime when other spin 
relaxation mechanisms due, e.g., to hyperfine and spin- 
orbit interactions are suppressed. We note that Eq. (5) 
is independent of Fl, whereas Eq. (6) is independent of 
Fi,;, so an asymmetric coupling of the DQD to the leads 
will play a role in determining the relative scales of the 
two contributions. 

At high temperature (/cbT > |w|) we have F{uj,T) ~ 
kgT, a regime that has been explored previously [13, 15]. 
In this work, we focus on the low-temperature regime 
(fcsT < |a;|), where F{uj,T) ~ uj9{uj), giving rates that 
vanish linearly for small energy separation, with signif- 
icant consequences (allowing, e.g., for the initialization 
of a pure spin state - see below). In Eq. (6), wc have 
assumed |A — e| 3> {wapl and have neglected resonant 
cotunneling contributions [20], which are exponentially 
suppressed for A/2 > fc^T. Corrections due to spin ex- 
change with the right lead are smaller in A/U' <C 1, 
where U' is the interdot charging energy. Additionally, 
we have considered the resonant tunneling regime [21] 
ie<t). 

We have numerically solved for the stationary current 
using the rates given in Eqs. (4)-(6) and have plotted the 
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result vs. B = b/g^B in Fig. 2. There is a sharp cutoff 
in the leakage current at large b {\b\ > \Es_\), which can 
be understood directly from Fig. 1(d). When the lowest- 
energy triplet state (T_ for 6 > 0) is below the lowest- 
energy singlet {S-), current will be blocked as soon as r_ 
is populated, since the transition from r_ to 5_ vanishes 
as Wt_^S- oc liJt_.s_0{ujt_,S-)- Thus, at low temper- 
ature a pure spin state can be prepared (|T+) or |T_) 
depending on the sign of b). We note that this prepara- 
tion can be achieved even in the presence of a large bias 
A/z > ksT . This is a nontrivial result, since a large bias 
will generally drive the DQD out of equilibrium, result- 
ing in a stationary state that is not determined by ther- 
mal equilibrium with the leads [22]. Moreover, using this 
method a pure spin state can be dynamically prepared 
on a time scale Tprep ~ (A/r^)^ ~ 0.1 /is (using pa- 
rameter values from the caption of Fig. 2) without the 
need to wait for slow spin relaxation processes. 

In Fig. 2, we show cuts at e = describing the evo- 
lution of I{B) as the field gradient 5B = Sb/gfiB is in- 
creased from zero (see the auxiliary material [23] for the 
dependence on e). For sufficiently large SB^ a dip ap- 
pears near B = 0. Similar zero-field dips have been seen 
experimentally in several DQD systems and have been 
attributed to effects due to hyperfinc [7] or spin-orbit 
[11, 24] coupling. In the present context, this zero- field 
dip can be understood from Fig. 1(d), without additional 
spin relaxation mechanisms. When Sb (or Tn) is large. 
To has a fast direct escape path by virtue of Eq. (5), so 
only the r+ and T_ states can block current. At 6 = 0, 
all triplets are degenerate, resulting in a vanishing inelas- 
tic cotunneling rate at low temperature (Wtj-^To — 0); 
transport can only occur if T± escapes via However, 
for a small nonvanishing Zceman splitting 6 > 0, we have 
Wtj^—,To oc 5 7^ 0, allowing an additional escape route for 
T-f. This results in an initial rise in current for small 6, 
which eventually must fall to zero when b ~ V2t, where 
r_ goes below S-. In the limit where Wto^u ^ Wa^p, 
we find a simple expression for the stationary leakage 
current: 



A 



(V2t-\b\){V2t + ?,[ 
V2t + \b\ 



^-eiV2t-\b\), (7) 



where we have taken e = 0, and T = 0. The current 
reaches a maximum at femax = ±V2t(2/A/3— 1) ~ ±0.22i. 
Eq. (7) is shown as a dash-dotted line in Fig. 2. We note 
that the hmit Wtq^o ^ Wa~^p required for Eq. (7) can 
also be achieved for much smaller 6b when Tr ^ Tl- A 
sufficiently large electron temperature will wash out the 
zero-field dip, but provided 6max ^ ksT, this feature will 
still be visible (see the inset of Fig. 2). Reaching this 
regime for T ~ 100 mK should be possible in nanowire 
[24] or nanotube DQDs [25] where t > 100 fj,eV is com- 
mon. 

We now turn to the transient (time-dependent) current 
that fiows between blocking events. We consider the in- 



stant after an electron has tunneled from the DQD to the 
right lead. With spin-independent tunneling rates, this 
leaves the dot in an equal mixture of the states |t) and 
II), setting the initial condition: Pa{0) = 1/2, Pa{0) = 0. 
The transient current into the right lead is then given by 
lB.{t) ~ ^^a^a-^<jPa{t)- The average number of elec- 
trons m that passes through the DQD, given a charge 
collection (measurement) time Tm is 



m{TM) = - 



dTlnir) 



(8) 



In Fig. 3 we plot m{TM) found by integration of Eq. 
(3) for a range of parameters when the stationary cur- 
rent is zero (i.e., ksT = e = 0, 5 > \/2t). The ac- 
cumulated charge shows a series of plateaux at time 
scales determined by the three types of rates given in 
Eqs. (4), (5), and (6). To better understand these 
plateaux, we consider the long-time saturation value 
TO = limTjvf-*oo "^(Ja/), which has been measured exper- 
imentally [2] and can be evaluated directly. 

We assume a probability Pb for the DQD to be in 
a blocking state each time an additional electron tun- 
nels onto the DQD. The probability of exactly n elec- 
trons passing through the DQD before current is blocked 
is then P„ = (1 — Pb)"^Pb, from which we find m = 
J2n = (1 ^ Pb)/Pb- In the simplest case, there may 
be Nb blocking levels out of N total, giving Pb = Nb/N . 
Assuming the two spin-polarized triplets \T±) are perfect 
blocking states, we set Pb = 1/2 since 2 out of 4 (1, 1)- 
states block perfectly, giving to = 1 (the expected value 
in Ref. [2]; the measured value was m ~ 1.5). How- 
ever, in the presence of some decay mechanism, tran- 
sitions between the various two-electron states (see the 
inset of Fig. 3) no longer allow for a clear definition 
of "blocking" levels. Nevertheless, we can still deter- 
mine Pb from the sum of probabilities for each inde- 
pendent path leading to a blocking state (T_ for the 
case shown in Fig. 3): Pb = Pa~*T- + Pa~*To^T- + 
Pa^T+^To^T- = (1 + P + 9P)/4, where Pa~.b... indi- 
cates the probability for a transition from state A to B, 
etc., and where the branching ratios are given (for & > 0) 
by: p = Wto^t_/(Wto^t_ +Wto-*s_ +Ea ^To^-y) and 
q = WT+^Tal{WT+^To+WT+^S-+WT+^s+)- Inserting 
this result gives: 



TO 



p-pq 



1 + p + pq 



(9) 



The average effective charge transported between block- 
ing events e* = (to -|- l)e is non-integral in general, rang- 
ing from e* = |e to e* = 4e. Eq. (9) allows for a 
precision measurement of slow spin-relaxation processes 
characterized by p and q, independent of the microscopic 
mechanism [26]. For concreteness, we consider the effects 
of cotunneling and field-assisted sequential tunneling be- 
low. 
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FIG. 3: (Color online) Average number m{TM) of electrons 
passing through the DQD within measurement time Tm- All 
values are as in Fig. 2 with the addition of b = l.Olv'St, 
T = 0, and Vr = IOFl = 10 /xel/. m{TM) is shown for 
(5B = OmT (solid line), (5B = 3mT (dash-dotted line), and 
5B = 100 mT (dashed line). The predicted saturation points 
for 3 of 4 levels blocking (m = 1/3), 2 of 4 levels blocking 
(m = 1) and 1 of 4 levels blocking (m = 3) are shown with 
dotted lines. The decay cascade (inset) defines the branching 
ratios p and q. 



In the limit of zero detuning (e = 0), we find q = 1/2, 
independent of b and t, leaving 



m 



m 



V2t + 2\b\ 



ASbV 

TLt J 



(10) 



r. 



Thus, at e = 0, m can be tuned from m = 3/5 to m = 3 
by varying 6, Sb, T^^fj, and t. Eq. (10) correctly predicts 
the saturation values at m=0.96, 1.4, and 3.0 for the 
solid, dash-dotted, and dashed lines, respectively, in Fig. 
3. 

We have analyzed the effects of inelastic cotunneling 
and a magnetic field gradient on the PSB. We find and ex- 
plain a zero-field dip in the stationary current, which may 
help to explain recent experimental results [7, 11, 24]. We 
have shown that a pure spin state can be dynamically ini- 
tialized, even at large bias, which is an important step on 
the way to full control over the quantum states of elec- 
tron spins. We have offered a possible explanation for an 
anomalously large value of the effective charge passing 
through the DQD found in experiments [2], which is im- 
portant for single-spin resonance studies. Our expression 
for this effective charge can be used to probe slow spin 
relaxation processes in the DQD to help understand the 
underlying physical mechanisms. A fractional effective 
charge e* in transport is often taken as evidence of ex- 
otic electronic states [27]. Here, we have shown that e* 
can reach universal fractional values in a simple system, 
without many-body correlations. 
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AUXILIARY MATERIAL FOR: "STATIONARY 
AND TRANSIENT LEAKAGE CURRENT IN 
THE PAULI SPIN BLOCKADE" 

Identifying spin decay 

In this supplement we show how to extract the 
magnetic-field dependence of microscopic relaxation pro- 
cesses at low field from the observed number of electrons 
that pass through the double dot between blocking events 
in the Pauli spin blockade regime. This measurement can 
be used to distinguish between spin-orbit, hyperfine-, and 
cotunneling-mediated spin relaxation mechanisms at low 
magnetic field, where other methods for single-spin de- 
tection fail. [SI] 

We recall the definitions for the branching ratios p and 
q (given before Eq. (9) of the main text). 



P 

q 



Wto-.T- 



Wt+^Tq 



(Sl) 
(S2) 



These expressions can be simplified under certain exper- 
imental conditions. In particular, we consider the case 



(a) 




(b) 




when there is no magnetic field gradient present, Sb = 0, 
and a sufficiently weak Zeeman splitting, so that tran- 
sition rates from triplets to triplets (split by the Zee- 
man energy) are smaller than those from triplets to sin- 
glets (split by exchange) {Wt„,^t^, ^ Wt„^s±)- With 
Sb = 0, the field-assisted sequential-tunneling rates van- 
ish: 



Wt 



0, 



and for a generic spin-flip Hamiltonian [S2], 



(S3) 



(S4) 



The rate Wto-»S_ is independent of Zeeman splitting, 
b = gfiB, since the Tq — splitting is independent of 
the global field. In contrast, Wtt oc b^ for small mag- 
netic field, depending on the spin-flip mechanism (as we 
have shown in the main text, rj ~ 1 for cotunneling- 
mediated spin-flips at low temperature, and previous 
work has shown 77 = 5 for spin-orbit-mediated spin flips 
with phonon emission [S3], and 77 = 3 for hyperfine- 
mcdiated spin fiips with phonon emission [S4]). For suf- 
ficiently small Zeeman splitting b, the triplets become 
degenerate resulting in vanishing rates, allowing us to 
approximate 



Wtt 

Wt+-,s± 



■P 



Wtt 

Wto-*S- 

Wtt 



Wt„~,s 
Wtt 



<C 1 
< 1 



< 1 



(S5) 
(S6) 
(S7) 



In this regime, we approximate m by its leading-order 
form in the small ratio Wtt /Ws±'- 



m = 3-2 



3-p~pq ^ 
1 + p + pq 

Wtt 



2p + 0{p^) + 0{pq), (S8) 

2\ 



O 



Wtt 
Wt^ 



■s 



(S9) 



Thus, by measuring m{b) = 3 — 76'', it is possible to ex- 
tract the relevant spin-flip mechanism: rj ~ 5,3, or 1 for 
spin-orbit interaction with phonon emission [S3], hyper- 
fine interaction with phonon emission [S4], or cotunnel- 
ing, respectively. 



PROCESSES LEADING TO DOT EXCITATION 





FIG. SI: ja) and are double dot energy eigenstates and 
AE = Ea — E/3 > 0. a) Available cotunneling processes b) 
Neglected cotunneling processes 



In our analysis we have neglected processes that can 
lead to excitation of the double dot at finite bias (see, 
e.g.. Fig. Sl(b) for an example). These processes can 
lead to nonvanishing stationary populations of excited 
dot states, which would correspond to "initialization er- 
rors" in the scheme we have proposed. However, we find 



6 




e (meV) 



FIG. S2: Magnetic field and detuning dependence of leakage 
current in the spin-blockade regime. We have taken T = 
40 mK, /ii = /ifl = 10 fj.eV , t = 100 fieV , A = 1 meV, SB = 
200 mT, and g = 2.0. A density plot shows a suppression in 
the current at _B = and a sharp cutoff when \gfiBB\ > \Es_ \ 
(solid white lines). The current runs from 1 = (dark blue) 
to / = 1.7pA (dark red). 



tion of the order (-0) which is negligible in our chosen 
regime. 

DETUNING AND FIELD DEPENDENCE OF 
CURRENT 

In Fig. S2, we show a density plot of the stationary 
current as a function of magnetic field and detuning for 
typical experimental parameters. The current shows a 
suppression at _B = corresponding to the zero-field dip 
in Fig. 2 of the main text. Solid white lines are drawn 
to indicate when the lowest-energy triplet state becomes 
the ground state (i.e., when \giJ,BB\ > |£'5_|). Fig. 2 of 
the main text corresponds to a cut along e = of Fig. 
S2, indicated here with a white dashed line. 



that the leading excitation processes are suppressed by 
the small parameter 

where U is the energy cost for double occupancy of one 
of the dots. Specifically, we have neglected virtual transi- 
tions to (1, 0)-charge states. Here, W^^"^ is the transition 
rate from |a) to |/?) due to process a{b) in Fig SI. The 
purity of the initial state will be reduced by a correc- 
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